An improved cosmological bound on the thermal axion mass 
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Relic thermal axions could play the role of an extra hot dark matter component in cosmological 
structure formation theories. By combining the most recent observational data we improve previous 
cosmological bounds on the axion mass m a in the so-called hadronic axion window. We obtain a 
limit on the axion mass m a < 0.42 eV at the 95% c.l. (m a < 0.72 eV at the 99% c.L). A novel aspect 
of the analysis presented here is the inclusion of massive neutrinos and how they may affect the 
bound on the axion mass. If neutrino masses belong to an inverted hierarchy scheme, for example, 
the above constraint is improved to m a < 0.38 eV at the 95% c.l. (m a < 0.67 eV at the 99% c.L). 
Future data from experiments as CAST will provide a direct test of the cosmological bound. 
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INTRODUCTION 

Recent Cosmic Microwave Background and Large 
Scale Structure surveys such as WMAP and SDSS have 
open the possibility of constraining fundamental physics 
with cosmology (see e.g. [H, Important upper lim- 
its on neutrino masses and energy densities, for example, 
have been obtained which are in some cases one order 
of magnitude better than the corresponding laboratory 
constraints (@, B 0, H) or competitive with big bang 
nucleosynthesis constraints (@j). 

The cosmological limits are model dependent and 
therefore rely on the assumption of a theoretical model of 
structure formation that, even if in agreement with cur- 
rent data, may need further key ingredients to explain 
mysteries and inconsistencies such as dark energy. More- 
over, for some datasets, the relevance of systematics is 
still matter of debate. 

However, future laboratory experiments will certainly 
test the cosmological results. The overlap of cosmological 
and laboratory limits will open a new window of investi- 
gation and may provide evidence for new physics and/or 
improve our knowledge of systematics. 

It is therefore timely to constrain fundamental physics 
with cosmology. In this paper we indeed move along one 
of those lines of investigation, providing new bounds on 
the thermal axion mass from cosmology. There are two 
possible ranges of axion masses (~ fxeV and ~ eV) and, 
in principle, both could provide either a dominant or a 
sub-dominant dark matter component. Here we focus on 
thermal axions with masses of ~ eV. For a recent revival 
of the cold dark matter scenario with axions of masses 
~ ^eV, see Ref. (zj. New constraints on the thermal 
axion mass and couplings have recently been presented 
by the CAST experiment, which searchs for axion- like 
particles from the Sun which couple to photons [8] . While 
the axion mass region probed by the CAST experiment 
is one order of magnitude lower than the cosmological 



bound presented here, an overlap of the two results is 
clearly around the corner. 

Let us remind the origin of the axions. Quantum 
Chromodynamics (QCD) respects CP symmetry, despite 
the existence of a natural, four dimensional, Lorentz 
and gauge invariant operator which badly violates CP. 
The former extra CP violating-term gives rise to phys- 
ical observables, namely, to a non- vanishing neutron 
dipole moment, d n . The existing tight bound \d n \ < 
3 x 1CP 26 ecm Q requires the CP term contribution to 
be very small. Why are CP violating effects so small in 
QCD? Why is CP not broken in QCD? This is known as 
the strong CP problem. The most convincing, and ele- 
gant, solution to the strong CP problem was provided by 
Peccei and Quinn by adding a new global U(1)pq 
symmetry. This symmetry is spontaneously broken at 
a large energy scale f a , generating a new spinless parti- 
cle, the axion, allowing for a dynamical restoration of the 
CP symmetry. Axions are the pseudo Nambu-Goldstonc 
bosons of the broken U(1)pq symmetry 11, 12| and may 
be copiously produced in the early universe, either ther- 
mally 13] or non-thermally [14j . providing a possible 
(sub)dominant (hot) dark matter candidate. The axion 
mass and couplings are inversely proportional to the ax- 
ion coupling constant f a 
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where R = 0.553 ± 0.043 is the up-to-down quark masses 
ratio [l5| and f n = 93 MeV is the pion decay con- 
stant. In principle, axions can interact with photons, 
electrons and hadrons. If axions couple to photons and 
electrons, the simplest bound comes from an energy loss 
argument. The axions produced in a star escape car- 
rying away ener gy, producing anomalous stellar observ- 
ables, see Refs. }16l 17 . [HI for a review. However, in 
practice, axion interactions are model dependent. Here 
we focus on hadronic axion models such as the KSVZ 
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model [19j, |20[, in which there is no tree level interac- 
tion between axions and lcptons and the axion-photon 
coupling could accidentally be negligibly small. Hannes- 
tad et al (2lj have recently found an upper limit on the 
hadronic axion mass m a < 1.05 eV (95% CL), which 
translates into f a > 5.7 x 10 6 GeV. In this letter, we re- 
inforce the former limit by means of an updated analysis, 
using a broad set of the most recent available cosmolog- 
ical data, and allowing for two possible hot dark matter 
components: neutrinos and axions. 



rameter H(T), becomes larger than their thermally aver- 
aged interaction rate. To compute the axion decoupling 
temperature To we follow the usual freeze out condition 



T{T D ) = H{T D ) . 
The axion interaction rate T is given by [2^ 
T = n^^nirij (a ljV ) , 



(4) 
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where n a = ((3/tt 2 )T 3 is the number density for axions 
in thermal equilibrium, and the sum extends to all pro- 
duction processes involving as initial states the particles 
i and j , which are in equilibrium at To . We will assume 
that the axion decay constant f a is sufficiently small to 
ensure that axions decouple from the thermal plasma af- 
ter the QCD transition epoch at T = T QC d ^ 200 MeV 
(fa < 4 x 10 7 GeV, i.e., m a > 0.14 eV). Consequently, 
we do not have to consider axion interactions with the 
quarks and gluons before the QCD phase transition and 
the dominant processes contributing to the thermally av- 
eraged cross section in Eq. |5]) will be 7r°7r ± — > a-K^ and 
7r + 7r~ — > air , see the interaction lagrangian, Eq. {2]). 
(■7T°TT + d^TT^ + 7T°7r _ 9 A1 7r + — 2n + tt~ d^n ) ,(2)We follow here the computation carried out by Chang 

and Choi [22| for the average rate ir + ir — > tt + a: 



THE HADRONIC AXION MODEL 

Among axion couplings with hadrons, those of interest 
for us are the axion-nucleon couplings C a N, responsible 
for the processes N + N <-> N + N + a and N + tt <-> 
N + a, and the axion-pion couplings £ a7r , responsible for 
a + 7r <-> 7r + tt. In practice, nucleons are so rare in the 
early universe respect to pions, that only the axion-pion 
interaction will be relevant for thcrmalization purposes. 
The lagrangian reads Q 



where 
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where 



is the axion-pion coupling constant [22| . The most strin- 
gent limits on the axion-nucleon coupling in hadronic ax- 
ion models, g a N — Czvmiv/ f a , are those coming from 
SN 1987A neutrino data. If axions couple to nucleons 
strongly, the supernova cooling process is modified, dis- 
torting both the measured neutrino flux and the dura- 
tion time of the neutrino burst emitted. The limit in the 
axion-nucleon coupling g a N, assuming that the model- 
dependent parameter CV — 0(1), translates into an ax- 
ion decay constant f a < few xl0~ 6 GeV [23| . Even if 
axion emission does not affect the SN cooling, if g a N 
is strong enough, the axion flux may excite 16 O nu- 
clei in water Cherenkov detectors. The absence of a 
large signal from radiative decays of excited 16 O* nu- 
clei in the Kamiokande experiment provides a lower limit 
f a > 3 x 10 5 GeV [24j. In summary, hadronic axions with 
the decay constant f a around 10 6 GeV, i.e. m a ~ eV, can 
escape from all astrophysical and laboratory constraints 
known so far, suggesting an ideal hot dark matter candi- 
date, within the mixed hot dark matter scenario 251. 



AXION DECOUPLING 

Axions will remain in thermal equilibrium until the 
expansion rate of the universe, given by the Hubble pa- 
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Here f(y) = l/(e y — 1) denotes the pion distribution 
function, Xl = \pi\/T, Vl = Ek/T (i = 1,2), s = 2{m\ + 
T 2 (yiy2 — X1X2U))), and we assume a common mass for 
the charged and neutral pions, = 138 MeV. 

The RHS in Eq. (gj) contains the Hubble expansion 
rate, related to the e nerg y density of the universe via the 
Friedmann equation [141 ]: 



H(T) 



— m— 
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where M p i is the Planck mass. We have computed, for 
temperatures T in the range 1 MeV < T < 200 MeV, 
i.e. between BBN and the QCD phase transition eras, 
the number of relativistic degrees of freedom g*(T), ac- 
cording to Ref. 14]. We neglect the axion contribution 



to g+ for simplicity. After resolving the freeze out equa- 
tion Eq. (@J|, we obtain the axion decoupling temperature 
To versus the axion mass m a (or, equivalently, versus 
the axion decay constant f a ). From the axion decou- 
pling temperature, we can compute the current axion 



number density, related to the present photon density 
ri 7 = 410.5 ± 0.5 cm" 3 E3 via 
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(9) 



where g+s refers to the number of entropic degrees of 
freedom. Before electron-positron annihilation at tem- 
peratures ~ eV, the number of entropic degrees of free- 
dom is g*s = g*, since all relativistic particles are at 
the same temperature. At the current temperature, 
g*s (T ) = 3.91 Q. 



COSMOLOGICAL CONSTRAINTS 

As now common practice in the literature we derive 
our constraints by analyzing Monte Carlo Markov Chain 
of cosmological models. For this purpose we use a modi- 
fied version of the publicly available Cosmo-MCMC pack- 
age cosmomc [26[ with a convergence diagnostics done 
through the Gelman and Rubin statistic. We sample the 
following eight-dimensional set of cosmological parame- 
ters, adopting flat priors on them: the baryon and Cold 
Dark Matter densities, ujf, — Slbh 2 and lu c — fl c h 2 , the ra- 
tio of the sound horizon to the angular diameter distance 
at decoupling, 6 S , the scalar spectral index rag, the over- 
all normalization of the spectrum A at k = 0.05 Mpc -1 , 
the optical depth to reionization, r, the energy density 
in massive neutrinos 



n v h 2 = 



92.5 eV ' 

and the energy density in the thermal axions: 



(10) 
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1.054 • 10 4 eV cm" 3 131 eV \ gi , s (T D ) 

(11) 

where we have used Eq. ©. For instance, for the 
hadronic axion upper mass bound quoted in Ref. [2l| . 
i.e. m a ~ 1.05 eV, the axion decouples at T D ~ 64 MeV, 
at which g+s{T D ) - 15.24 and Q a h 2 ~ 0.0053. 

We consider a combination of cosmological data which 
includes the three-year WMAP data [l|, t he small-scale 
CMB measurements of CBlJpTj, VSA j28|, ACBAR [H 
and BOOMERANG-2k2 [§(|7ln addition to the CMB 
data, we include the constraints on the real-space power 
spectrum of red luminous giant (LRG) galaxies from the 
fourth data release of the SLOAN galaxy redshift survey 
(SDSS) [1H and 2dF [H, and the Supernovae Legacy 
Survey data from [33| . Finally we include a prior on the 
Hubble par ameter from the Hubble Space Telescope Key 
project [34| and the BBN prior in form of a Gaussian 
prior on Qbh 2 (see e.g. Q). We refer to this dataset as 
Conservative in the rest of the paper. 



Conservative 


1.4 / 2.0 eV 0.55 / 0.9 eV 


Conservative+LYA 
+Sm„ > 0.05eV 
+Sm„ > O.leV 


0.42 / 0.72 eV 0.20 / 0.37 eV 
0.41 / 0.71eV 0.22 / 0.38eV 
0.38 / 0.67eV 0.25 / 0.44eV 



Conservative+LYA+BAO 0.35 / 0.64 eV 0.18 / 0.31 eV 

TABLE I: This figure shows the 95%/99.9% upper confidence 
limits on the marginalised posterior probabilities for axion 
and neutrino masses. See text for discussion. 



In the second dataset we include constraints on the 
small scale linear power spectrum coming from Lyman-a 
analysis of SDSS quasar spectra (3?], [HI . 

The main results of our analysis are reported in the 
Table HI As can see, without assuming any prior on the 
neutrino mass, the mass of the thermal axion is found 
to be m a < 0.42 eV and the sum of the three active 
massive neutrinos E m v < 0.20 eV, both at the 95% c.L, 
i.e Vt a h 2 < 0.0014 and tt u h 2 < 0.0018. Therefore, the 
neutrino-axion (hot) dark matter contribution represents 
a small fraction (< 2.5%) of the total CDM. Excluding 
from the analysis the constraints from BAO and Lyman- 
a cosmological datasets the former limits translate into 
rn a < 1.4 eV and E m u < 0.55 eV. The inclusion of the 
Lyman-a data has an enormous impact on the analysis. 
In the same table we also consider the effect of adding 
Baryonic Acoustic Oscillations (BAO) data detected in 
the Luminous Red Galaxies (LRG) sample of the SDSS 
[3r^ to the data. Strictly speaking this is a statistically 
incorrect procedure as the correlations with SDSS LRG 
power spectrum are not well understood, but it gives 
the idea of the improvements that can be achieved by 
including BAO constraints. 

In Fig.[T] where we present marginalized constraints on 
the E m v ~ m a plane. There is a clear anti-correlation 
between the constraints on the thermal axion mass and 
the mass of the three active neutrinos. In other words, 
the cosmological data allow only for a very specific quan- 
tity of hot dark matter: if one increases the active neu- 
trino mass, more hot dark matter is present in the model 
and the axion mass has to be smaller in order to fit the 
observations. 

Figure [2] depicts the 95% CL axion mass limits in 
the ma-gaj-y (axion-to-photon coupling) plane. The lim- 
its should be within the region allowed by the KSVZ 
model. We have considered two possible scenarios, ac- 
cordingly to neutrino oscillation data: normal hierarchy 

(E m - 



> 



'13 



> 



0.05 eV) and inverted hierarchy 



(E m ^ ~ 2 • ^/|Amf 3 | > 0.1 eV), as well as the massless 
neutrino case. The 95% c.L constraints that we obtain 
for the axion mass within the three possible scenarios 
mentioned above are m a < 0.34,0.31 and 0.34 eV, re- 
spectively, including both BAO and Lyman-a datasets. 
We found no significant difference between the normal 
hierarchy and the massless neutrino scenarios. If fu- 



FIG. 1: Likelihood contour plot in the ^2m v — m a plane 
showing the 68% and 95% c.l. from the conservative dataset 
(left panel) and from the complete dataset (right panel) . Note 
different axes. 



ture cosmological data [40( or direct terrestrial searches 
for neutrino masses, as the ones which will be carried 



out by the KATRIN experiment [4l|, improve the cur- 
rent limits on ^m„, one could obtain automatically a 
rather robust, independent, albeit indirect limit on the 
axion mass m a . We depict in Fig. [2] the current 95% 
c.l. CAST limit for comparison The CAST exper- 
iment has been upgraded and in the near future it will 
explore QCD axions, that is, a range of axion masses up 
to about 1 eV. Cosmology-independent future limits on 
the axion mass are therefore extremely important, since 
they could provide a test of the cosmological constraint 
and be translated into a limit of the universe's hot dark 
matter fraction in the form of massive neutrinos. 




(eV) 



FIG. 2: 95% CL limits on the axion mass obtained in the con- 
servative and full analysis (shaded regions), assuming three 
possible values of the sum of the neutrino masses in the m a - 
ga-y-y plane. From right to left the region represent the ex- 
clusion limits assuming a prior Yl fn v > 0, m v > 0.05 eV 
(N.H.) and > 0.1 eV (I.H). As a comparison, we show 

the recent results from the CAST experiment (blue contour) 
as well as the theoretical KSVZ parameter region (within the 
green lines), following Fig. 8 from Ref. and the CAST 
prospects (blue dashed line) [39l |. 
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CONCLUSIONS 

We have presented an improved limit on the hadronic 
axion mass by combining the most recent available cos- 
mological data. A novel content of this analysis is the 
addition of a hot dark matter component in the form of 
massive neutrinos. Interestingly, we have noticed an anti- 
correlation between the thermal axion mass and the mass 
of the three active neutrinos ^2 m v ■ This anti-correlation 
is due to the suppression induced on the small scale power 
spectra by both the relic axion and the massive neutrino 
free-streaming species. A larger (smaller) axion mass 
content can be traded by a smaller (larger) massive neu- 
trino content. R the complete cosmological dataset is 
used, we find m a < 0.35 eV and m v < 0.17 eV at the 
95% c.l., implying that the fraction of (hot) dark mat- 
ter in the form of massive thermal axions and neutrinos 
is only a few percent (< 2.5%) of the total CDM con- 
tent. The former limits get modified if priors on the neu- 
trino or axion masses are imposed. Future cosmological 
and/or terrestrial searches for neutrino (axion) masses 
could therefore be translated into an improved and inde- 
pendent axion (neutrino) mass limit. 
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